Surface Flows of Granular Mixtures: III. Canonical Model 
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We present the generalization of the minimal model for surface flows of granular mixtures, 
proposed by Boutreux and de Gennes [J. Phys. I France 6, 1295 (1996)]. The minimal model was 
valid for grains differing only in their surface properties. The present model also takes into account 
differences in the size of the grains. We apply the model to study segregation in two-dimensional 
silos of mixtures of grains differing in size and/or surface properties. When the difference in size is 
small, the model predicts that a continuous segregation appears in the static phase during the filling 
of a silo. When the difference in size is wide, we take into account the segregation of the grains in 
the rolling phase, and the model predicts complete segregation and stratification in agreement with 
experimental observations. 

83.70.Fn - Granular solids 

83.10.Pp - Particle dynamics 

47.55.Kf - Multiphase and particle-laden flows 



SHORTENED TITLE: Surface flows of granular mixtures 
I. INTRODUCTION 

Segregation of mixtures of grains is commonly observed 
in granular materials that are poured, vibrated, or ro- 
tated [ [pL| — P^9|] ] - This phenomenon has great importance 
in industrial processes. The simplest way to observe seg- 
regation is to pour a mixture of grains of different sizes 
onto a heap; one obtains a heap with the large grains near 
the bottom and the small grains at the top [ pQ-Eq]]. 

The mixture can also be poured in a two-dimensional 
silo made up of two vertical plates separated by a gap of 
approximately 5 mm (a granular Hele-Shaw cell), as stud- 
ied recently by Makse et al. [ p3 Different forms of 
segregation are observed when the grains differ both in 
their size and their surface properties (shape, roughness, 
stickiness). When the grains do not differ much in size, 
the larger or the smoother grains stop preferentially at 
the bottom of the slope, and the smaller or the rougher 
at the top of the slope. This segregation is limited since 
grains of both species remain present everywhere; this 
phenomenon is called continuous segregation. 

When the grains have a wide difference in size, seg- 
regation is stronger and appears in two different ways. 
When the large grains are smoother than the small ones, 
a complete segregation is observed, where all large grains 
stop near the bottom of the slope, and all small grains 
stop near the top. When the large grains are rougher 
than the small grains, a spectacular stratification is ob- 
tained, where the grains deposit in alternating layers of 
different species, parallel to the sand-pile surface. 

Theoretical studies of surface flows of grains were trig- 
gered by the works of Bouchaud, Cates, R. Prakash, and 
Edwards (BCRE) [ and Mehta and collaborators [ 

[B2fl]. BCRE proposed a set of coupled equations to de- 



scribe granular flows in one s pec ies sand- piles. Boutreux 
and de Gennes (BdG) [ ||,||] generalized the BCRE 
equations in order to describe granular mixtures com- 
posed of two species. They proposed a theoretical formal- 
ism, and a minimal model describing the case of grains 
with different surface properties but equal size. In this 
case, BdG found a power law behavior of the concentra- 
tions which explains the continuous segregation. Since 
the minimal model does not take into account the size 
difference between the grains, Boutreux has treated the 
important case where the grains differ only in size, in a 
second article [ |3^,|3(|] of the series started by BdG [ 

Hi- 

In the present article, the third and last one of the 
series, we treat the general case of the canonical model, 
by describing different examples where the grains differ 
in size and/or surface properties. We first assume that 
the two species do not differ much in size. We then con- 
sider that there is no segregation inside the rolling phase, 
which is homogeneous in the vertical direction. In this 
case, we show that the canonical model predicts continu- 
ous segregation (in agreement with experiments) with a 
power law behavior of the concentrations. 

In the last part of the article, we consider the case 
in which the two species have a wide difference in size. 
We argue that in this case segregation probably occurs 
directly inside the rolling phase; due to percolation of 
the small grains, the small grains fall through the gaps 
between the large ones. We show that by taking this phe- 
nomenon into account and relaxing the assumption that 
the rolling phase is homogeneous, complete segregation 
occurs when the small grains are rougher than the large 
grains, and stratification occurs when the large grains are 
rougher than the small grains, in agreement with exper- 
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iments. In these cases, our results are consistent with 
theoretical studies recently published by Makse, Cizeau, 
and Stanley (MCS) [ f3?]-|3j|], who proposed a modified 
version of the BdG equations. MCS proposed a model 
valid for a mixture of grains with very different size. They 
explained the complete segregation phenomenon by pre- 
dicting an exponential behavior of the concentrations, 
and successfully reproduced the mechanism leading to 
stratification as observed in the experiments. 

The present article is organized as follows. In Sec. ||, 
we review the theoretical formalism developed in [ 
po| Q]. In Sec. Ill, we present our model of binary colli- 



sions between one rolling grain and one grain at rest, and 
we obtain the canonical model for surface flows of granu- 
lar mixtures. Section IV describes the application of the 



model to all possible cases, when the two granular species 
differ in size and/or in surface properties. Then in Sec. 
we model the steady state filling of a two-dimensional silo 
by a mixture of grains, and we discuss the predicted seg- 
regation profiles according to the different composition of 
the mixture. Finally in Sec. VI, we consider the case in 
which the grains have a wide difference in size, in order 
to describe stratification and complete segregation. 



rolling phase 




of the static phase (see Fig. |l|). We call 9{x,t) the local 
slope of the interface, and h(x, t) the height of the static 
phase. For notational convenience we do not consider the 
difference between the angle 9 and its tangent 



6(x,t) 



dh 

dx 



(1) 



We call 4> a {x, t) the volume fractions of the two species 
of grains in the static phase just below the interface (here 
the index a denotes the grain species, and is equal to T' 
or '2'). We have (f>\ + 4>i — 1. We assume that both 
species have a small difference in size, so that there is 
no segregation inside the rolling phase, i.e., the rolling 
phase is homogeneous in the vertical direction (this as- 
sumption will be released in Sec. VI). We call R(x,t) the 
total thickness of the rolling phase. We also consider two 
"equivalent thicknesses" for the two species in the rolling 
phase R a (x 1 t) (i.e., the total thickness multiplied by the 
local volume fraction of the a grains in the rolling phase 
at position x). The total thickness of the rolling phase is 
then equal to 



R(x,t) = Ri{x,t) + R 2 (x,t). 



(2) 



The equations describing surface flow of grains take 
into account the fact that grains in the rolling phase move 
downwards due to their weight, and collisions between 
the rolling grains and the static grains induce exchanges 
between the two phases. The equation that describes the 
exchange of grains between the two phases is 



-(R 



1 coll 
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(3) 



FIG. 1. Diagram showing the variables used to describe the 
granular flow of mixtures in a two-dimensional silo. 



exchange between the a grains in the rolling phase with 
the static phase. Equation (||) can also be written for the 
single species 



-R a 



coll ■ 



(4) 



II. THEORETICAL FORMALISM 

The study of two-dimensional surface flows of granular 
materials h as p rogressed significantly with the works of 
BCRE [ HQ] and Mehta et al. [ §§]]. BCRE proposed 
a set of variables and a set of coupled equations to de- 
scribe two-dimensional granular flows for a pure species of 
grains. Recently, BdG [ J34|] generalized the BCRE for- 
malism by considering a mixture of two granular species. 
In the present paper, we describe this formalism for two 
species. Following Bouchaud et at, we assume that there 
is a sharp distinction between a static phase, where grains 
at rest belong to the pile, and a thin rolling phase where 
grains are not part of the pile but roll downwards on top 



The evolution equation for each species in the rolling 
phase, taking into account the downhill convection of 
grains due to gravity, is 

Ra = V ~Q^~ + Ralcolh (5) 

where v is the convection speed of the rolling grains, 
which may depend on the slope 9(x, t) and on the grain 
species. However, we assume that the two species are 
mixed inside the rolling phase; the convection speed v 
must be identical for the two species. Moreover, in prac- 
tice variations in the angle 9 are small enough, so that 
the speed v can then be taken as a constant. 
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FIG. 2. Diagram showing the different type of collisions 
between a rolling grain of type 1, and a static grains of type 
1 or 2. oi: auto-amplification, X\\ cross-amplification, 6i: 
auto-capture, z%: cross- capture. 



III. CANONICAL MODEL 

We propose a microscopic model of the grain collisions, 
that will allow us to calculate the exchange term R a \coii 
as a function of 9, R a , and <p a . We calculate, in first order 
approximation, the interaction term R a \ C oii by consider- 
ing the binary collisions between one rolling grain and 
one grain at rest on top of the static phase. Collisions 
between three grains (or more) are less probable, and will 
be neglected. We consider four types of collisions for a 
rolling grain of type 1 (Fig. |J): 

• Auto-amplification: another type 1 grain starts 
to roll. This collision contributes with a term 
ai(9)(f>iRi to Ri\coii- This term is proportional to 
i?i because all type 1 rolling grains interact with 
the static phase due to the fact that the rolling- 
phase is thin. It is also proportional to the con- 
centration of type 1 grains in the static phase, <fii . 
The function ai(9) is called the collision function 
for the auto-amplification of type 1 grains. This 
function is positive, has dimensions of a frequency, 
and depends a priori on the angle 9. The collision 
functions may also depend on the type of species in 
the static phase that are in contact with the grain 
about to roll; the rougher the neighboring grains, 
the less chance for the static grain to start rolling, 
and the smaller the value of a\. However, this 
neighbouring effect is expected to be small com- 
pared to the variations of a\ with respect to 9. In 
order to simplify, we will neglect the neighbouring 
effect, and will consider that ai only depends on 
9. The same approximation will be used for all the 
other collision functions. 



• Cross-amplification: a type 2 grain starts to move. 
This collision contributes with a term xi(9)<p2Ri to 

R2\coll- 

• Auto-capture: the type 1 rolling grain is captured 
after a collision with a type 1 static grain. This pro- 
cess contributes with a term —bi(9)<j)iRi to R\\ C oil- 

• Cross-capture: the type 1 rolling grain is captured 
by a type 2 static grain. This process contributes 
with a term — zi(9)4>2Ri to i?i| co ;;. This cross- 
capture interaction was not taken into account in 
the minimal model [ [M]. It plays an important 
role when the grains have different sizes, as will 
appear later. 

Collisions where the rolling grain stops and simultane- 
ously a static grain starts to roll are a limiting case be- 
tween cross-amplification and cross-capture. The proba- 
bility that these two collisions happen simultaneously is 
low, and therefore we neglect this process. 

When the colliding rolling grain belongs to the type 2 
species, four similar binary collisions occur. We call the 
corresponding collision functions 0,2(9), X2(9), &2($), and 
Z2{9), which are all positive. Since increasing the slope 
of the pile favors rolling, the amplification functions a a 
and x a are increasing functions of 9. Conversely, since 
decreasing the slope favors capture, the capture functions 
b a and z a are decreasing functions of 9. 

It is now possible to write the expressions of the ex- 
change term R a \ C oii in a matrix form. We define the 
collision matrix M by [ fl3~4| ] 

\R2\coi1) \R2j y ' 

The previous microscopic description yields 
M=( ^- b ^' Zlh , u& A - (?) 

V El 02 (0-2 - 02)02 ~ Z 2 <pl J 

The elements M Q ^ are determined by the concentrations 
(f> a (x,t), and by the local angle 9(x,t) via the collision 
functions. 

Next we calculate the expressions of the collision func- 
tions by doing a linear approximation around the angles 
of repose. In order to do so, let us consider the quantity 
Ei defined as the total exchange from the static phase 
to the rolling phase due to collisions originated by type 
1 rolling grains — similar definition can be provided for 
the type 2 rolling grains, i?2- A type 1 rolling grain 
can interact — via auto-amplification or auto-capture — 
with other type 1 static grains giving rise to a contribu- 
tion (ai — bi)4>iRi to Ei, or it can interact — via cross- 
amplification or cross-capture — with type 2 static grains 
giving rise to a contribution (x± — zi)4>2Ri to E±. Then 
Ei is given by 

Ei = [(ai - &i)0i + (xi - Zi)<t> 2 ]Ri 

= (Mn + M 2 i)Ri. (8) 
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Notice that E% contributes to Ri\ co u via M%x, and also 
to R 2 \ co a via M 2 \. 

In a model of a single species made of type 1 grains, 
the exchange term would be equal to E\ — \a\(G) — 
b\(ff)\R\(x, t), and the angle for which a\ = b\ would 
correspond to the situation where there is no exchange 
of grains between the static and the rolling phases. We 
call this angle the angle of repose of the pure type 1 
grains, which we denote On. In general, the repose an- 
gle 6 aa of the pure a species is the angle at which the 
auto-amplification and auto-capture functions intersect 



(9) 



Moreover, we will call the cross-angle of repose a p the 
angle for which the cross-amplification function equals 
the cross-capture function 



X a (0afj) = z a {Q a p). 



(10) 



The angle 0\ 2 for the type 1 grains is defined by X\ (#12) = 
zi(0\ 2 ), and corresponds to the angle of repose of a type 
1 rolling grains moving on top of a surface of type 2 static 
grains [ |3^]]. 

When the surface properties of the two species are not 
very different, the two pure angles of repose do not dif- 
fer very much. Moreover, when the size of the grains do 
not differ much, the cross-angle of repose are also close. 
In practice, the angle 9 remains close to the angles of re- 
pose of the species and we can linearize the eight collision 
functions with respect to 9 to get a simple expression for 



Ei = hi{0-0 n )fa +~/ 2 (6 -0i 



]Rx, 



(11) 



where 



71 = deai - dgbi, 72 = de%i - de^i- (12) 

In order to consider the simplest case, we assume that all 
the derivatives of the collision functions have the same 
order of magnitude, i.e. 

deai ~ dex\ ~ —dgbi ~ —dezi, (13) 

so that 7! ~ 72 = 7, and therefore 

E 1 = 7 [9-0 1 (fa)]R 1 , (14) 

where the angle Oi(fa) is given by 

(15) 



n\<P2) 



712 - 0\\)(p2. 



The constant 7 has the dimensions of frequency, and has 
typical value 7 ~ 25 s _1 [ ^7j]. Dimensional analysis 
show that the order of magnitude of 7 is given by 



7 ~ v/d, 



(16) 



where d is the typical size of the grain. The constant 7 
represents the frequency of interaction between a rolling 



grain and the static phase. The larger the value of 7, the 
more frequent the exchange between phases. 

Equation fli~4| ) shows that 0\(fa) is a cross-over an- 
gle; Ei describes capture of rolling grains (E\ < 0) when 
9 < 0i, and amplification of grains (E\ > 0) when 9 > 9\. 
The angle 0\(fa) for a mixture of grains plays the role 
of the constant angle of repose 0\\ for a pure species of 
grains; the angle Q\(fa) is called the generalized angle of 
repose for the type 1 species. When no type 2 grains are 
present on the static phase (fa = 0), 0i(fa) (given by Eq. 
(|l5|)) is equal to the angle of repose 6*11 of the pure type 

1 species. When fa increases, the type 1 rolling grains 
interact more frequently with type 2 static grains, and 
the generalized angle of repose 0\(fa) changes. 

The gen eralized angle of repose was introduced by 
MCS [ |37j]], who have shown that it has a key role in 
explaining the segregation as well as the stratification of 
granular mixtures. MCS defined the generalized angle 
of repose a (4>p) of species a as a linear function of the 
volume fraction <j>p and proposed the expression ( |l5| ) . In 
the present paper we show that this expression can be 
derived, by using the binary collision model and a linear 
development approximation. When it is possible, nota- 
tions of reference [ |37|] are used in the present paper. 

Making similar approximations, we get a simplified ex- 
pression for the exchange E 2 due to collisions with type 

2 rolling grains 



where 



E 2 = ~f[9 ~ 9 2 (fa)]R 2 , 



l) 



(17) 



(18) 



is the generalized angle of repose for the type 2 species. 

The generalized angles of repose quantify the degree 
of interaction of the rolling species with the pile. If the 
species a have a smaller generalized angle of repose for 
any value of fa, then, when R\ = R 2 , the a rolling grains 
are less captured (or amplify more grains in the bulk), 
than the other species. In order to quantify this behav- 
ior, we define the difference 



12 



l). 



(19) 



In the general case, ipi 2 is a function of fa. However in 
order to simplify, in the following we will consider par- 
ticular situations where ^12 is a constant that does not 
depend on fa (i.e. 9 Y2 - 6>n = 22 - 9 2 \). 

It is now possible to write a simpler expression for the 
collision matrix (^). Previous calculations for E a yield 



M 



l[0-0i{fa)]-x 1 
xi(0)fa 



X 2 



l[0 - 02{fa)] - X 2 \U)(pi 

(20) 



where the cross-amplification functions can be written in 
the following way 
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zi(0) = MO - 0u) + x , x 2 (6)=x 1 (e)-Ax 12 . (21) 



Here Xq and Aa; 12 are two constants of the model. We 
call Eqs. ( p0| ) and ( fjl| ) the canonical form of the collision 
matrix. This canonical form is general, and valid for any 
mixture of two granular species. In practice, the grains 
differ in size and in surface properties. Let us now see 
how the canonical form of the collision matrix conveys 
these differences. 
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FIG. 3. Case A: grains differing only in their sizes (small 
type s and large type I grains), (a) Diagram of the collision 
functions, (b) Diagram of the generalized angles of repose as 
a function of the volume fractions <j>i- The angles of repose 
0ss and Bu for the two pure species are equal. Due to the size 
difference, the cross-angles of repose 9i B and 8„i are different. 



IV. POSSIBLE DIFFERENCES BETWEEN THE 
GRANULAR SPECIES 



A. Grains differing only in their size 

Let us consider the case where we mix a species of 
small grains, denoted the V species, with a species of 
large grains, denoted the 'V species, both species having 
similar densities and the same shape or roughness (in this 
case, indices '1' and '2' are replaced by indices V and 'I'). 
The size difference between the two species allows us to 
compare the collision functions. A grain more easily sets 
a small grain into motion than the reverse, hence 



x s (9) < xi{0). 



(22) 



A small grain is more easily captured on a surface of large 
grains than the reverse, then 



zi(6) < z s (9). 



(23) 



Equations 



and 



defined by Eq. 



|23|) imply that the two cross-angles 
(nQ), satisfy 



Hs < Vsl- 



(24) 



The size difference between the two species does not al- 
low to compare neither the functions ai{9) and a s (9), nor 
the functions bi(9) and b s (9). 

However considering that the 'I' and V species have 
the same surface properties, we can obtain other inequal- 
ities. In first order approximation, the probability that 
an V static grain is set into motion by an 1 s 1 rolling 
grain is equal to the probability that an '/' static grain is 
set into motion when it is collided by an 'V rolling grain. 
The two probabilities are approximately equal because in 
both collisions the two interacting grains have the same 
weight, and because the two interacting surfaces are iden- 
tical. Hence we have 



a; 



(25) 



Similarly, the probability that a rolling grain is captured 
by a static grain belonging to the same species does not 
depend on the grain species being ' V or then 

b s ~ h = b. (26) 

Equations ( |25| ) and (^6|) imply that the two pure angles 
9 aa (defined by Eq. (g)) satisfy 

9 u = 9 SS . (27) 

Thus, the present model shows that two granular 
species with identical surface properties have also identi- 
cal angles of repose, as observed in experiments. Finally, 
the size difference between the two species involved in 
collisions yields 



x s < a < xi, zi < b < z s 



(28) 



We consider the simplest case by using the set of simplest 
relations consistent with Eqs. (|28|) 



a = (x s +xi)/2, b=(z s + zi)/2. 



(29) 
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Equations fl2q ) imply the following relations for the an- 
gles of repose, as proposed by MCS [ J37J] 

9 ls < 6 U = 9 SS < 9 sl . (30) 



into motion. The larger the size difference between the 
two species, the larger the value of ip s i and the larger the 
strength of the exchange processes. In our case, ip s i is 
a small constant so that the linear development of the 
collision functions remains valid. 
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FIG. 4. Case B: grains differing only in their surface prop- 
erties (rough type r and smooth type m grains), (a) Collision 
functions, (b) Generalized angles of repose: due to the sur- 
face differences the two angles of repose 8 3S and 9u of the pure 
species are now different, but the cross-angles of repose are 
the same due to the equal size of the grains. 

The collision functions and the generalized angles of 
repose 9i(4> s ) and # s (</>;) are represented in Figs. ||a 
and ||b. Note that due to our assumptions (p^), ip s i = 
6 s (4>i) — Qi{4>s) is a constant independent of cf> a . The an- 
gles of repose 8 SS and 9u of both pure species are equal 
due to their equal shapes, but due to the size difference 
between the particles we have ip s i > for any value of 
(j) a ; the I rolling species amplifies more easily the rolling 
phase than the s rolling species, because large grains are 
less easily captured, and more easily set another grain 



B. Grains differing only in their surfaces properties 

Let us now consider the case where we mix species 
of rough grains, denoted the V species, with species of 
smooth grains, denoted the 'm' species. The rougher the 
surface of a grain species, the more efficient a collision 
between two grains belonging to this species. Hence we 
have 



a m (8) < a r (6), b m {9) < b r (9). 



(31) 



Experiments show that granular species with the 
smoother surface have also the smaller angle of repose: 
9mm < 9 rr [ p7|]. This inequality is equivalent to the 
condition 



dm ^ b r 



(32) 



Thus the last inequality is an experimental constraint for 
our model. 

The different surface properties between the two 
species does not allow to compare neither the functions 
x m and x r , nor the functions z m and z r . However, if 
we also assume that the V and 'm! species have the 
same size, both cross-interactions correspond to collisions 
where two grains of the same mass interact via a rough 
surface in contact with a smooth surface. In first order 
approximation, the probability of a cross- interaction does 
not change if the two grains involved in the collision are 
switched, hence 



Equations (|33|) imply 



y rm — u mr- 



(33) 



(34) 



Finally, the surface differences between the two species 
involved in collisions yield 

a m < x < a ri b m < z < b r . (35) 

We take the simplest relations consistent with Eqs. ( |35| ) 

x = (a m + a r )/2, z=(b m + b r )/2. (36) 



Equations ( J35| ) imply the following relations for the an- 
gles 



< 



< 



(37) 



The collision functions and the generalized angles of re- 
pose 9 r ((j)m) and 9 m (<f> r ) are shown in Figs. [|a and [|b. 



G 



Here again ip rm is a small constant due to our assump- 
tions (|3q ) and to the validity of the linear approxima- 
tions. Since m rolling grains are less easily captured, 
they amplify more easily the rolling phase than the r 
rolling grains. 
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FIG. 5. Case C: mixture of small rough grains (type 1) and 
large smooth grains (type 2). (a) Collision functions: when 
the grains differ both in size and in shape, the eight functions 
are distinct, (b) Generalized angles of repose. 



C. Mixture of small rough grains and large smooth 
grains 

We denote "type 1" grains the small rough grains, and 
"type 2" the large smooth grains (here we have 1 = s = r, 
and 2 = I = m). The size difference between the two 
species still yields Eqs. (|22|) and ((2^), and the cross-angles 
of repose verify 621 < Moreover, the surface differ- 
ences still imply Eqs. ([H]), and the angles of repose ver- 
ify 022 < If the size difference is very small, effects 
due to the surface differences are stronger, and we get 



#22 < #21 < #12 < 0ii- This case is close to the situation 
where grains differ only in their surface properties. 

However, when the size difference is more important, 
we obtain a case close to the situation where the particles 
differ only in their sizes. Then, we have (as postulated 
in [ ©]) 

021 < 022 < 011 < 012, (38) 



and 



x\ < a 2 < ai < X2, Z2 < 62 < &i < X\. (39) 



The collision functions and the generalized angles of re- 
pose for this case are represented in Figs. ||a and ||b. 
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FIG. 6. Case D: mixture of small smooth grains (type 1), 
and large rough grains (type 2) when the size difference is 
small, (a) Collision functions, (b) Generalized angles of re- 
pose. 
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D. Mixture of small smooth grains and large rough 
grains 



- w = 7 [0 - 0i(0a)]iZi + 7 [0 - e 2 ^))R 2 . (45) 



We denote "type 1" grains the small smooth grains, 
and "type 2" the large rough grains (1 = s = m and 
2 = I = r). As explained previously, we have for the 
different angles of repose 



'11 < #22, 



721 < °12- 



(40) 



We consider here the case when the difference in size 
is small, and shape segregation effect is more important 
than size segregation effect. Then in this situation, we 
have (see Fig. gb) 

011 < 021 < 012 < 022- (41) 

As shown in Fig. ^a, the collision functions satisfy 

a\ < x\ < x 2 < a 2 , b\ < z 2 < Z\ < b 2 . (42) 



pouring point 




L x 



V. SEGREGATION IN THE FILLING OF A SILO 

A. General equations 

In order to describe the segregation between the two 
species predicted by the present model, we study the case 
of the steady state filling of a two-dimensional silo, which 
is a simple geometry and has been the focus of experimen- 
tal studies by different authors [ J|^-|29|]]. Inside a two- 
dimensional cell made of two vertical plates separated by 
approximately 0.5 cm, a mixture of two different species 
is poured with a constant in-going flux. Let us call the 
horizontal axis the x axis. The cell is located between 
< x < L, and the grains are poured at the point of 
injection at x = L (Fig. 0). We focus on the steady state 
filling; as the mixture is poured in the cell, the surface of 
the pile rises uniformly without deforming, at a constant 
speed w 



dh 



(43) 



We review some results found for this situation with 
the minimal model [ In the steady state filling, 

we have R a = 0. Then equations (||) and (||), and the 
boundary condition R(x = 0) = 0, imply that the to- 
tal thickness R(x) of the rolling phase decreases linearly 
with respect to the distance from the pouring point 



R(x) 



w 



-x. 



(44) 



The total exchange between the two phases R\\ co u + 
R2\coii can be calculated as a function of the generalized 
angles of repose, using the collision matrix (ppf). Then 
Eq. (|) yields 



FIG. 7. Filling of a two-dimensional silo, in which a con- 
stant flux of grains is poured. In the steady state, the surface 
of the static grains rises uniformly, at the constant speed w. 

In this expression, the order of magnitude of the right 
hand side is ji/jR ~ ipwx/d, where 7 ~ v/d [Eq. (|6|)]. 
The ratio of the right hand side of ( |45| ) divided by its left 
hand side is ~ ?px/d. Thus, when x 3> d/ip, the left hand 
side can be neglected, and we get for the local slope 



Ri 
R 



Ra 
R 



(46) 



In particular, in the steady state, the value of lies be- 
tween the two generalized angles of repose 0i and 9 2 . By 
combining Eqs. (|j), (^0|) and (f46|), we obtain the follow- 
ing expressions of the volume fractions tj> a in the static 
phase 



'i = ( 1+ fia 



R2 



Ri 

x\R\ + x 2 R 2 J R ' 
Ri \ R2 



X1R1 + x 2 R 2 J R ' 



where 



£12 = 7^12 + Ax 12 



(47) 



(48) 



is a constant in our model. We notice that the minimal 
model of [ |Q] consists of taking the cross-amplification 
functions constant independent on the angle, x\{6) = 
x 2 (9) — m, and Axi 2 — 0. With these approximations 
the minimal model can be solved to obtain a closed form 
for the concentrations and the rolling grain profiles. 

Equations d47| ) show that segregation happens when 
rolling grains stop on the static phase; the volume frac- 
tion R a jR of the a grains in the rolling phase is differ- 
ent from the volume fraction <p a of the same species in 
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the static phase. Since the quantity R a /{x\Ri + X2R2) 
is always positive, if £12 > we obtain <f>i > Ri/R 
and 4>2 < R2/R', type 1 grains are captured more eas- 
ily than type 2 grains. Conversely, if £12 < we obtain 
4>\ < Ri/R and 02 > R2/R', type 2 grains stop more 
easily. This behavior induces segregation everywhere in 
the static phase; the larger | £12 | the stronger the segre- 
gation. 

For instance, let us consider the case where we pour 
an equal volume mixture, R\(x = L) = R2(x = L). If 
£12 > 0, we obtain at the top of the slope (j>i(x = L) > 
4>2{x — L). Type 1 rolling grains stop more easily, and 
the fraction of these grains decrease in the rolling phase 
(i?i < R2), while type 2 rolling grains will preferentially 
stop at the bottom of the slope (j>i (x — 0) < ^{x = 0). 
Thus, in this example, the static phase contains mostly 
type 1 grains in its upper part, and type 2 grains in its 
lower part. 

In the lower part of the pile (x <C L) 1 it is possible to 
quantify more precisely the segregation for any kind of 
poured mixture. By doing a linear development of the 
equations, we can look for a power law behavior valid at 
the bottom of the slope, as suggested by the solution of 
the minimal model [ jj4|]. If, for instance, we assume £12 
to be positive, we expect that Ri/R tends to zero. We 
write for x <C L 



Ri{x) 
R(x) 



Ax", 



(49) 



where A and 77 are two positive constants. Then <f>i can 
be calculated in first order approximation with Eqs. ( |47| ) 



0i (x) 



^ ) .!,-". 

X2 



(50) 



When we substitute these expressions for Ri/R and (f>i 
in the relation dRi/dx = w<p\/v (obtained from Eqs. (^) 
and @), we obtain an equation whose compatibility con- 
firms that the type 1 species follows a power-law behavior 
at the bottom of the slope. This equation also gives the 
value of the exponent rj 



V = 



6 



X2 ( 



7 2 -2 



(51) 



Note that the exponent rj depends on the granular 
species. This power law behavior shows that our model 
predicts continuous segregation; the concentrations of 
species vary slowly in the container, both types of grains 
remaining present everywhere. This continuous segre- 
gation corresponds to what is observed in experiments, 
when the two granular species do not have a wide differ- 
ence in size [ 29 ] . 




FIG. 8. Case B: grains differing only in their surface prop- 
erties. Concentration profiles (calculated numerically) during 
the steady state filling of a silo. The segregation clearly ap- 
pears at x = L, where R m = R r but <f> m < <f> r . 



B. Segregation for the different mixtures 



For the different mixtures of grains treated in Sec. IV, 
the model predicts continuous segregation during the 
steady state filling of a silo. In the case of grains dif- 
fering only in their surface properties, we have x r = x m 
and 



7Vv 



Aa 



> 0. 



(52) 



The sign of £ rm indicates that inside the static phase, 
the rough grains are found preferentially at the top of 
the silo [see Eq. fl47|)]. In order to quantify more pre- 
cisely the predictions of the model, we perform a nu- 
merical integration of the equations of motion for this 
case. Figure || shows the results calculated in the case 
of an equal volume mixture of rough grains and smooth 
grains, R r (x = L) = R m {x — L). For this simula- 
tion, we chose 9 rr = 49°, 9 mm = 43°, tfj rm = 3°, and 
x(6 = 6 rm ) = 0.3757. Figure || shows the volume frac- 
tions <p a (x) for < x < L in the steady state. Note that 
segregation clearly appears at x = L, where R m = R r 
but <j) m < 4> r . Figure ^ confirms that our model predicts 
a continuous segregation, and not a complete one; r 
does not fall suddenly at x = L/2, but slowly decreases 
as x decreases. Similar profiles have been recently ob- 
tained with numerical simulations, by using a granular 
media lattice gas model to study the filling process of a 
two-dimensional silo with inelastic particles differing in 
friction coefficients [ p0|]. 

When the grains differ only in their size, a similar 
continuous segregation is found, with the smaller grains 
located preferentially at the top of the pile. Thus our 
results indicate that the roughest grains or the small- 
est ones are segregated at the top of the pile, while the 







largest grains or the smoothest ones are segregated pref- 
erentially at the bottom. When we mix grains differing 
both in size and in surface properties, a competition be- 
tween size-segregation and shape-segregation appears. 




FIG. 9. Case C: small rough grains (type 1) and large 
smooth grains (type 2). Concentration profiles at the sur- 
face calculated numerically in the steady state. The dotted 
line corresponds to the fit to the analytical expression ([so|), 
valid at the bottom of the pile. 



For the mixture treated in Sec. [V C| of small rough 
grains and large smooth grains, the model predicts con- 
tinuous segregation with the rough and small grains 
found preferentially at the top of the pile. This is the case 
of strongest segregation found with the present model, 
since both size segregation and shape segregation act si- 
multaneously to segregate the small rough grains at the 
top of the pile and the large smooth grains at the bottom. 
Figure [| shows the volume fractions <pa (%) calculated nu- 
merically when an equal volume mixture of small rough 
(type 1) grains and large smooth (type 2) grains is poured 
in the silo. We set x = OA-f, Ax 12 = -0-137, L = 100, 
and use the collision functions shown in Fig. [| Figure ^ 
shows a continuous segregation which is stronger than in 
the case o f grai ns differing only in their surface proper- 
ties (case 
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Fig. |3|). At the lower end of the slope 
(i « L), we have a complete purification of both species 
due to segregation; Ri(x)/R and 4>i(x) tend to zero. In 
that region, we can fit the numerical solution with the 
analytical expression Eq. ( |50|) with r/ = 0.29. 

The last case treated in Sec. ^ corresponds to a mix- 
ture of small smooth grains and large rough grains. The 
size difference is taken to be small, to allow for the lin- 
ear development of the collision functions. Thus shape 
segregation is found to dominate over size segregation, so 
that the smooth grains are found at the bottom. How- 
ever, when the effect of size segregation is comparable to 
the effect of shape segregation, the competition between 
size and shape segregation gives rise to a different phe- 
nomenon. Indeed, experiments show that we obtain ei- 



ther stratification when the large grains are the roughest, 
or a complete segregation when the large grains are the 
smoothest [ [|6| ||)]. The theoretical model proposed in 
the present article can explain both the complete segre- 
gation and the stratification. These phenomena are due 
to the type of segregation that appears directly inside the 
rolling phase, as explained in the next section. 



VI. COMPLETE SEGREGATION AND 
STRATIFICATION 

So far we have treated the cases when the difference 
between surface properties and size of the grains is not 
too wide, hence we could perform linear developments 
of the collision functions around the region of interest. 
This situation gives rise to the continuous segregation 
patterns in all the cases studied in the previous sections. 
When the difference between the grains size is important, 
stronger segregation effects are expected. According to 
experiments [ p6|-p9|], this occurs when p > 1.5, where p 
is the ratio of the size of the large grains divided by the 
size of the small ones. The size segregation is due to seg- 
regation at the shear surface between the rolling and the 
static phase [ |p7fll , or it may already happen inside the 
rolling phase [ |2 7| , |30|l 1 ; small rolling grains tend to fall 
downward through the gaps in between the large grains, 
so that they form a sub-layer of small rolling grains un- 
derneath the large rolling grains. This phenomenon is 
called "kinematic sieving", "free-surface segregation" or 
"percolation" [ |2^^,^]]. The large rolling grains are 
not in contact with the bulk, and are captured after the 
small grains. The percolation effect takes place inside 
the rolling phase if its thickness is larger than (approxi- 
mately) two or three grain diameters. However, even for 
a thin flow, strong size segregation occurs in the shear 
surface between the rolling and static phase if the size 
difference between the grains is wide enough. In both 
cases, the large rolling grains are captured only after the 
small ones, and the collision functions may take the forms 
proposed in [ j37|] (see Fig. 4 of [ @]). 

Here we use a suitable modification of the interaction 
term, in order to take into account the segregation in- 
side the rolling phase within the present formalism. We 
replace R 2 (x,t) in the definition of the interaction term 
by R 2 (x, 1) exp[Ai?i(x, t)/R(x, t)]. The exponential fac- 
tor multiplying R 2 (x,t) mimics the fact that the inter- 
action of large rolling grains is screened by the presence 
of small grains, so that large rolling grains R 2 interact 
with the grains at the surface of the static phase only 
when Ri(x,t) <§C R(x,t)/X. The dimensionless parame- 
ter A > measures the degree of percolation [ We 
numerically simulate this model using a mixture of grains 
differing both in size and shape. When the large grains 
are smoother than the small grains, we obtain a complete 
segregation of the mixture (see Fig. |l0|a). The transition 
zone between the two species at rest has a size which goes 
approximately from a few mm to a few cm. When the 
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large grains are rougher than the small grains, Fig. [To| b 
shows that we find stratification of the static phase. This 
stratification results from the competition between size- 
segregation and shape-segregation. In conclusion, when 
our model includes the percolation effect, we are able to 
reproduce the experimental observations obtained with 
particles of very different sizes [ p6|-p9|]. Our results 
are also consistent with the ones found by Makse et al. [ 
§7| H], t0 whom 

we refer for further details. 




x L 




FIG. 10. (a) Concentration profiles obtained numerically 
showing the complete segregation of a mixture of small rough 
grains (type 1) and large smooth grains (type 2) when the 
percolation effect is present, (b) Stratification of a mixture 
of large rough grains (dark) and small smooth grains (grey) , 
when the percolation effect is included in the model. 



VII. DISCUSSION 

We have proposed an analytical model (called canon- 
ical model) that explains the continuous segregation ob- 
served during the filling of a silo, when the grains have 
a small difference in size. Our predictions could be pre- 
cisely tested by experimental measurements. When the 
two granular species have a wide difference in size, we 
include, in the canonical model, the percolation effect 
that appears directly in the rolling phase. The model is 
then able to reproduce the complete segregation and the 
stratification observed in experiments; it is also consis- 
tent with results recently published by Makse et al. who 
described grains with a wide difference in size. 



The canonical model incorporates both differences in 
size and in surface properties of the grains, and con- 
tains several coupling parameters between the two gran- 
ular species such as the cross-angles of repose, and the 
two constants xn and Ax±2- In comparison, the min- 
imal model [ |Q] is simpler than the canonical model 
since the minimal model contains only one coupling pa- 
rameter. However, as a consequence, the minimal model 
cannot take into account a size difference between the 
particles. 

In the present paper, we have applied the canonical 
model to only one practical situation: the steady state 
filling of a silo. This situation is interesting, since it 
has already been described experimentally by different 
authors. Moreover, the present model for surface flows 
of granular mixtures could be applied to many different 
situations. In particular, the segregation that appears 
during thick avalanches in rotating cylinders could be 
studied analytically; the results could be compared to 
the ones published recently for a pure granular species [ 
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